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Abstract 

We present a method for tabulating all cubic function fields over ¥q{t) 
whose discriminant D has either odd degree or even degree and the leading 
coefficient of —3D is a non-square in F,, up to a given bound B on deg(_D). 
Our method is based on a generalization of Belabas' method for tabulating 
cubic number fields. The main theoretical ingredient is a generalization 
of a theorem of Davenport and Heilbronn to cubic function fields, along 
with a reduction theory for binary cubic forms that provides an efficient 
way to compute equivalence classes of binary cubic forms. The algorithm 
requires O(B'^q^) field operations as B — >■ oo. The algorithm, examples 
and numerical data for q = 5,7, 11, 13 are included. 

1 Introduction and Motivation 

In 1997, Belabas [3] presented an algorithm for tabulating all non-isomorphic 
cubic number fields of discriminant D with \D\ < X for any X > 0. In the 
above context, tabulation means that all non-isomorphic fields with discrimi- 
nant \D\ < X are hsted or written to a file, by listing the minimal polyno- 
mial for each respective field. The results make use of the reduction theory 



for binary cubic forms with integral coefficients. A theorem of Davenport and 
Heilbronn |15) states that there is a discriminant-preserving bijection between 
Q-isomorphism classes of cubic number fields of discriminant D and a certain 
explicitly characterizable set lA of equivalence classes of primitive irreducible 
integral binary cubic forms of the same discriminant D. Using this one-to-one 
correspondence, one can enumerate all cubic number fields of discriminant D 
with \D\ < X hy computing the unique reduced representative f{x,y) of every 
equivalence class in U of discriminant D with \D\ < X. The corresponding 
field is then obtained by simply adjoining a root of the irreducible cubic f{x, 1) 
to Q. Belabas' algorithm is essentially linear in X, and performs quite well in 
practice. 

In this paper, we give an extension of the above approach to function fields. 
That is, we present a method for tabulating all Fg(i)-isomorphism classes of 
cubic function fields over a fixed finite field up to a given upper bound on 
the degree of the discriminant, using the theory of binary cubic forms with coef- 
ficients in Fq[t], where Fg is a finite field with char(Fg) ^ 2, 3. The discriminant 
D must also satisfy certain technical conditions, which we detail below. 

This paper is a substantial expansion and extension of the material in 
and corresponds to Chapters 4 and 6 of the first author's Ph.D thesis, prepared 
under the supervision of the last two authors. The present paper includes a 
more detailed description of the algorithm in the unusual discriminant case 
and improved bounds on the coefficients of a reduced form compared to those 
appearing in [26j . Another somewhat complementary approach to tabulating 
cubic function fields different from the one taken in this paper is described 
in [20], and relies on computations of quadratic ideals to construct all cubic 
function fields of a given fixed discriminant. The approach in [20] is limited to 
finite fields Fg where q = —\ (mod 3), and unlike our algorithm, considers only 
fixed discriminants and not all discriminants D with deg(Z3) bounded above. 

Our main tool is the function field analogue of the Davenport-Heilbronn 
theorem [15] mentioned above, which generalizes to Dedekind domains (see 
Taniguchi [3T]). As in the case of integral forms, we also make use of the 
association of any binary cubic form / of discriminant D over F^ \t\ to its Hes- 
sian Hf which is a binary quadratic form over Vq\t\ of discriminant —3-D. Under 
certain conditions on the discriminant, this association can be exploited to de- 
velop a reduction theory for binary cubic forms over Fg[t] that is analogous to 
the reduction theory for integral binary cubic forms. Suppose that deg{D) is 
odd, or that deg(Z3) is even and the leading coefficient of —31? is a non-square 
in F*. We will establish that under these conditions, the equivalence class of / 
contains a unique reduced form, i.e. a binary cubic form that satisfies certain 
normalization conditions and has a reduced Hessian. Thus, equivalence classes 
of binary cubic forms can be efficiently identified via their unique representa- 
tives. The case where deg(-D) is odd is analogous to the case of definite binary 
quadratic forms, but the other case has no number field analogue. 

Our tabulation method proceeds analogously to the number field scenario. 
The function field analogue of the Davenport-Heilbronn theorem states that 
there is a discriminant-preserving bijection between Fq(t)-isomorphism classes 
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of cubic function fields of discriminant D GVq [t] and a certain set U of primitive 
irreducible binary cubic forms over Vq[t\ of discriminant D. Hence, in order to 
list all Fq(t)-isomorpliism classes of cubic function fields up to an upper bound 
B on deg(-D), it suffices to enumerate the unique reduced representatives of all 
equivalence classes of binary cubic forms of discriminant D for all D g ^q[t] 
with deg(il') < B. Bounds on the coefficients of such a reduced form show that 
there are only finitely many candidates for a fixed discriminant. These bounds 
can then be employed in nested loops over the coefficients to test whether each 
form found lies in U. The coefficient bounds obtained for function fields are 
different from those used by Belabas for number fields, due to the fact that the 
degree valuation is non- Archimedean. In fact, we obtain far simpler and more 
elegant bounds than those in the number field case. 

This paper is organized as follows. Section [2] begins with some background 
material on algebraic function fields. The reduction theory for imaginary and 
unusual binary quadratic forms and binary cubic forms over [t] is developed in 
Sections [3] and |4j respectively. The derivation of the bounds on the coefficients 
of a reduced binary cubic form appears in Section [5] The Davenport-Heilbronn 
Theorem for cubic function fields is presented in Sectionl6] We detail the tabu- 
lation algorithm as well as numerical results in Section [7j Finally, we conclude 
with some open problems and future research directions in Section [8j 

2 Preliminaries 

For a general introduction to algebraic function fields, we refer the reader to 
Rosen [M] or Stichtenoth [30] • Let be a finite field of characteristic at least 5, 
and set F* = Fq\{0}. Denote by Fg[t] and Vq{t) the ring of polynomials and the 
field of rational functions in the variable t over Fg, respectively. For any non- 
zero H £ Fq[t] of degree n = deg{H), we let \H\ — q" — q^^^s^^), and denote by 
sgn{H) the leading coefficient of H. For H — 0, we set \H\ — 0. This absolute 
value extends in the obvious way to F, (t) . Note that in contrast to the absolute 
value on the rational numbers, the absolute value on ¥q{t) is non- Archimedean. 

An algebraic function field is a finite extension K of Fq(t); its degree is the 
field extension degree n = [K : ¥q{t)]. It is always possible to write a function 
field a,s K = ¥q{t,y) where F{t,y) = and F{Y) is a monic polynomial of 
degree n in F with coefficients in F, [t] that is irreducible over F^ (t) . We assume 
that ¥q is the full constant field of K, i.e. F{Y) is absolutely irreducible. 

A homogeneous polynomial in two variables of degree 2, with coefficients in 
F, [t] , of the form 

f{x,y) = Ax^ + Bxy + C?/^ 

is called a binary quadratic form over ¥q[t]. We abbreviate the form as / = 
{A,B,C). Similarly, a homogeneous polynomial in two variables of degree 3, 
with coefficients in F^ [t] , of the form 

f{x, y) — ax^ + bx^y + cxy^ + dy^ 
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is called a binary cubic form over Vg[t]. We abbreviate the form f{x,y) by 
/ = (a, b, c, d). The discriminant Z) of a binary quadratic form / — (A, B, C) is 
D{f) — — AAC. In a similar vein, the discriminant of a binary cubic form 
/ = (a, b, c, d) is D = D{f) = ISabcd + b^c^ - Aac^ - Ab^d - 21a^d'^. We wih 
assume that all forms are primitive, irreducible over Fg[i], and have distinct 
roots and thus non-zero discriminant. 

Let -D be a polynomial in [t] . Then D is said to be imaginary if D has odd 
degree, unusual if D has even degree and sgn(Z3) is a non-square in F*, and real 
if D has even degree and sgn(D) is a square in F*. Correspondingly, a binary 
quadratic form is said to be imaginary, unusual or real according to whether 
its discriminant is imaginary, unusual or real. These terms have their origins 
in quadratic number fields. If D is imaginary (resp. real), then the quadratic 
function field Fq(t, \/£)) shows many similarities to an imaginary (resp. real) 
quadratic number field, such as the splitting of the infinite place of Vq{t) and 
the unit group structure. For D unusual, there is no number field analogue to 
the function field Fg(i, \/Z))- The terminology "unusual" is due to Enge [TS] . 
and we adopt this terminology for quadratic fields and binary quadratic forms 
in this paper. 

Integral binary forms have a rich history going back to Lagrange and Gauss, 
and many important applications (see Buchmann and VoUmer |7j and Buell 
[5]). Their reduction theory was developed for positive definite forms first, as 
this case is the most straightforward. Recall (from Buell [8j, Chapters 1 and 2, 
for example) that an integral binary quadratic form {A^B^C) is definite if its 
discriminant is negative. In this case, both A and C have the same sign. One 
then further specializes to positive definite forms; these are forms with negative 
discriminant and A > (and hence C > 0). In other words, one considers 
the element / = (A, _B, C) in the associate class of definite forms with A > Q. 
Correspondingly, if / = [A^ B, C) is a binary quadratic form with imaginary 
or unusual discriminant, then we say that / is positive definite if sgn(A) is a 
square in F*, and negative definite otherwise. 

We summarize a few useful notions and results pertaining to binary quadratic 
and cubic forms over ¥q[t] below. These results are completely analogous to their 
well-known counterparts for integral binary quadratic forms, found in Chapter 
8 of Cohen 



{M : M is a 2 X 2 matrix with entries in ¥q[t] and det(Af) 7^ 0} 
is denoted by AIat2{¥q[t]). As usual, the set 



Let / be a binary quadratic or cubic form and M ~ \" s J ^ Mat2(¥q[t]). 
The action of M on / is given by f oM — fiax + Py, 'yx + Sy). Using this action, 
we give the definition of equivalence, which differs slightly from |26j in that the 



The set 



GL2{Vq[t]) = {M : M e Mat2{¥q[t]) with det(M) e FJ} 



denotes the general linear group of degree 2. 
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extra multiplication of / by a unit is removed. This was done to simplify some 
of the subsequent work. Two binary forms / and g over ¥q[t\ are said to be 
equivalent if 

f{ax + /3y, jx + Sy) = g{x, y) 

for some a,/?, 7, 5 £ '¥q[i\ with aS — fi^ S F*, i.e., g — f o M for some M in 
GL2{¥q[t]). We immediately obtain equations for the coefficients of equivalent 
binary quadratic forms, namely, if / = {A, B, C) and g — {A' , B' , C) — f o M, 
with M as above, then 

A' = Aa"^ + Ba-f + C7^ (1) 
B' = 2AaP + B{a6 + Pj) +2CjS, (2) 
C" = A/S^ + B(5S + C6^. (3) 

Analogous formulas for binary cubic forms are also easy to obtain, and are 
omitted. 

We note that it indeed possible for a positive definite unusual form to be 
equivalent to a negative definite form. For instance, if | A| = |C| with sgn(A) = 1 
(so / is positive definite) and g = 1 (mod 4) (so —h/A is a non-square), then 
swapping A and C yields an equivalent negative definite form. 

Proposition 2.1. Let M = (^"5) € Mat2{¥q[t]) and f a binary quadratic or 
cubic form over ¥q [t] . Then the following hold: 

1. D{f o M) — (aS — /37)^ • D{f), if f is a binary quadratic form. 

2. D{f o M) — (aS — (3^)^ ■ D{f), if f is a binary cubic form. 

3. If M is non-singular, then f o M is irreducible over Fg [t] if and only if f 
is irreducible. 

4. If M is non-singular, then f o AI is primitive if and only if f is primitive. 



By Proposition 2.1 up to an even power of det(M), equivalent binary 
quadratic and cubic forms have the same discriminant. In addition, primi- 
tivity and irreducibility are preserved by the action just defined, provided the 
transformation matrix is non-singular. 

We finish this section by introducing the Hessian of a binary cubic form, 
along with some of its properties which are easily verified by straightforward 
computation. 

Definition 2.2. Let f — (a,b,c,d) be a binary cubic form over ¥q[t]. The 
Hessian of f and the polynomials P, Q, R are given by 

Hfix,y) = --deti \^Px^ + Qxy + Ry\ 

\ dydx dydy / 

where P = b^ — 3ac, Q ^ be — 9ad and R = c^ — 3bd. 
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Note that Hf is a binary quadratic form over Vq[t]. 

Proposition 2.3. Let f = {a,b^c,d) be a binary cubic form over ¥q[t\ with 
Hessian Hf = {P, Q, R). Then the following properties are satisfied. 

1. For any M G Mat2{¥q[t]), we have HfoM = (det A/)^ ■ HfoM. 



2. D{Hf)^-3 -Dif). 

3 Reduction of Binary Quadratic Forms over 



In tlris section, we give the reduction theory for binary quadratic forms with 
polynomial coefficients, due to Artin 1 . A modified version of Artin's material 
is presented here, as he does not consider binary quadratic forms, but only 
their roots, which results in a simpler treatment. Furthermore, some of his 
presentation is streamlined in this paper, using more modern notation. 

The reduction theory for binary quadratic and cubic forms allows us to single 
out a unique representative in each equivalence class of forms. The theory also 
enables the efficient computation of this representative, and demonstrates that 
there are only finitely many such equivalence classes for any given non-zero 
discriminant D in the case of binary quadratic forms. The case where D is a 
real discriminant of a binary quadratic form is excluded from the paper. 

The following conventions will be adopted. As before, let ¥q be a finite field 
of characteristic at least 5. Fix a primitive root h of F*. We predefine the set 
S = {h' : < i < (q- 3)/2}, so that a e S" if and only if -a ^ S. As in Artin 
[T], the discriminant Z? of / is also endowed with the normalization sgn(Z3) = 1 
or sgn(Z)) = h, where 1 or /i is chosen depending on whether or not sgn(Z3) is 
a square in F*. We choose this normalization in order to avoid the possibility 
of forms being equivalent to each other while possessing different discriminants. 
By Proposition \2.1\ the discriminant of a form can only change by a square 
factor of a finite field element, so normalizing to a single square or non-square 
sign value accomplishes this task. 

We provide the reduction theory for binary quadratic forms over ¥q[t]. We 
will treat the imaginary and unusual scenarios in parallel where possible, but 
there are significant differences. The case of unusual binary quadratic forms 
differs from that of imaginary forms (see ^1 ) , as it has no number field analogue. 
Another crucial difference is as follows: the analogous definition of "reduced" 
does not lead to a unique representative in each equivalence class in the case 
that 1^1 = |C|, but instead to q -I- 1 equivalent forms, called partially reduced 
forms, defined below. To achieve uniqueness, a distinguished representative 
among these q + 1 equivalent partially reduced forms needs to be identified. 
Some of the results and proofs below apply to imaginary forms as well, and 
these occasions will be explicitly mentioned. 
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For the remainder of this section, let / = {A, B, C) be an imaginary or un- 
usual binary quadratic form with discriminant D. Recall that we only consider 
irreducible forms, hence A ^ 0. Furthermore, our earlier assumption that ¥g is 
the full constant field of the function field Fq(t, VJ)) implies that D ^¥q. 

Definition 3.1. An imaginary or unusual binary quadratic form f = {A, B, C) 
over ¥q [t] is partially reduced if it satisfies the following conditions: 

1. \B\ < \A\ < \C\; 

2. If\A\ < \C\, then sgn(A) G {l,h}, and if \A\ = \C\, then sgn(^) = 1; 

3. B ^ implies sgn{B) G S. 

Note that \B\ < \A\ < \C\ if and only if \B\ < \A\ < \D\^^'^, and that reduced 
unusual forms are always positive definite. The exponent in |Z3|^/^ — q'^°s{D)/2 
in the above lemma is a half integer in the case where / is an imaginary binary 
quadratic form, so property (1) above is in fact equivalent to \B\ < \A\ < \C\. 
However, in the case of unusual binary quadratic forms, the exponent in -\/|_D| = 
gdcg(£')/2 integer. Hence, equality |^| = |C| = may in fact occur in 

this case. 

Proposition 3.2. Every binary quadratic form is equivalent to a partially re- 
duced binary quadratic form f — (A, B, C) of the same discriminant. 

Proof. The procedure to find / satisfying property (1) above is completely anal- 
ogous to the one for integral binary cubic forms; see Algorithms 5.1-5.3, page 
87 of Buchniann and VoUmer [7], or [15]. Now, let / = {A,B,C) be a bi- 



nary quadratic form satisfying condition (1) of Definition 3.1 If \A\ < \C\ and 
sgn(^) ^ 1 01 h, then we replace /(x, y) by f{e~^x, ey), where = sgn(A) if / 
is positive definite and = sgii{A)h~^ if / is negative definite. In other words, 
replace / with / o M, where AI = ( '^^ o ) . This will yield a form satisfying 
\B\ < \A\ < \C\ and sgn(A) G {1, h}. 

If \A\ = \C\, then consider the norm map from ¥q{^/h) down to F^; this 
map is always surjective (see for example Theorem 2.28, p. 54 of 22 ). Thus, 
sgn(A)~^ G F* is the norm of some element a + eVh in ¥q{Vh) (a,e G ¥q). 
It is now easy to see from equations ([l])~([3]) that the transformation matrix 

(2sgn(A)e s^i{A)a^ of determinant 1 yiclds a binary quadratic form /' ~ {A' , B' ,C') 
with sgn(A') = 1 and \B'\ < \A'\ = \C'\. 

Finally, If i? 7^ 0, then sgn(i?) G is achieved by transforming / with 
J = ( J ) if necessary. □ 

If / = {A,B,C) is a partially reduced form with \A\ = \C\ and —1 is a 
non-square in F^, then —4/h — sgn(C)^^ is a square in F^, say —4/h = e^. 
Then {A,B,C) is equivalent to {e'^C, B,e~^A) via the matrix (° "^q^)- Hence, 
additional normalization conditions will be needed to obtain a unique represen- 
tative in each equivalence class in the case when \A\ = |C|. We now characterize 
when two partially reduced binary quadratic forms are equivalent. 
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Theorem 3.3. Let f = {A, B, C) and f = {A' , B' , C) he two partially reduced 
binary quadratic forms of the same discriminant. Then the following hold: 

1. If \A\ < \C\ and \A'\ < \C'\, then f and f are equivalent if and only if 

2. If \A\ < \C\ and \A'\ ~ |C"| then f and f are not equivalent. 

3. If \A\ = \C\ and \A'\ = \C'\ then f and f are equivalent with f'^foM 
if and only if 

a 13 
AuP/h ua 



M 



where a,/3 € F,, - {A/h)P^ = u := det(M) ^ ±1, and if B' ^ 0, 
then u is determined by the condition sgn(i3') e S. 

Proof. Suppose that /' = foM, with M = ( " ^ ) and m det(M) = a5~l3-f € 

¥*. Then u = ±1 by Proposition 2.1 since /, /' have the same discriminant. 

Since {A\B',C') = f{ax + (3y,-/x + Sy), equations hold. Since / = 

/' o M~^, we also obtain: 

uA = A'S^ - B'jS + C'j^, (4) 
uB = -2A'PS + B'{a5 + p-/)~2C'a-f, (5) 
uC = A'p^ B'aP + C'a^. (6) 

Forming appropriate products of ([T]), Q and ^ yields: 

4uAA' = {2Aa + B-f f - D-/^ , (7) 

4uCA' = {2C-1 + Baf - Da^ , (8) 

AuAC = {2A13 + B5f ~ D5^ , (9) 

AuCC = {2C6~Bpf~D(3^. (10) 

Suppose first that |^| < |C|, so \A\ < \D\^/^. Then by equation \A'\ < 
|i:)|i/2 and 

\{2Aa + Bjf - D-f^\ < \D\. 

Since the leading terms of the expressions {2Aa + B'-f)^ and Dj^ on the left hand 
side of the above inequality cannot cancel, this forces 7 = 0. Thus det(Af) = a6, 
and hence a,6 e ¥*. Now ^ implies \A\ = \A'\. Then \D\ = \AC\ = \A'C'\ 
precludes the possibility that \A'\ — \C'\, which proves part (2) of theorem. 
Furthermore, together with \B\, \B'\ < \A\ forces /3 = 0. Thus /3 = 7 = 0. 
Then equations (|l|-(|| yield A' = Aa^, B' = B{a5), C = 

It follows that det(M) — aS = ±1. Since both sgn(^) and sgn(A') are 
either squares (in which case they are both 1) or non-squares (in which case 
they are both /i), we thus obtain = 1. Hence a ~ ±.1. This yields four 
possibilities for the matrix M: ±/, ± J, where J = {q -i) ■ If -8 = = B' , then 
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f — f ■ If B ^ 0, then both sgn(i?), sgn(i3') <E S, which means that aS = 1. 
This leaves M = ±1, hence / = /', as desired, completing the proof of part (1). 

To prove part (3), suppose now that \A\ = \C\ = \A'\ = \C'\ = jZ?!^/^; in 
particular, D is an unusual discriminant. First, we deduce that a,P,j,5 G ¥g 
under this assertion. To see this, note that the absolute value of each of the 
left-hand sides of equations ([7])-([To]) above equals \D\. Note also that there 
cannot be cancellation of leading terms on the right-hand side of each of the 
above equations, since sgn(£>) is a non-square. Thus, the only way that ([7|-([l0| 
can hold is if a,/3,7, (5 G ¥q. Recall that sgn(^) = sgn(A') = 1, sgn(Z3) — h 
and thus sgn(C) — sgn(C") — —h/A. If we compare the coefficients of 
of both sides of Equations ([T])-([3]), we obtain 



af3 - 



-h 



1 : 

7^, 



h 



4 
+(5^ 



(11) 
(12) 

(13) 



If 7 = 0, then we deduce as before that M = ±1 or M — ±J, where 
J = ( Q _?i ) . These two cases can only occur when B = 0. 
If 7 ^ 0, then from equation (12), we obtain 

4 a/3 
h 7 



From equations (13) and (|11|), it follows that 

1 = 



16\ 



/l2 



r 



16/32 



Hence, 7^ {\Q/h'^)P^ and so 7 = ±{A/h)p. Write 7 = e(-4//i)/3 with e = ±1. 
Then 

= = —ea. 

fi 7 



Then by (11), 



Hence 



u — aS — — —e c? 



M 



-h 



e 7 = — e a 

4 ' ' 



a P 
{AuP)/h ua 



7 = -e. 



with a,f3 e ¥q, - {4/h)/3'^ = 1 and u = det(A/) = ±1 as desired. Any 
change in u between 1 and —1 clearly changes the sign of B' if i?' 7^ 0, so m is 
determined by sgn(i3') in this case. □ 
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We note that with the same notion of reducedness, Theorem 3.3 is not true 
for real binary quadratic forms, as one cannot deduce that a, f^,^,5 G ¥q in the 
same way as in the proof of Theorem |3.3| 



Part (3) of Theorem 3.3 now yields the following: 



Corollary 3.4. Any unusual partially reduced binary quadratic form f = (A, B, C) 
satisfying \A\ — \C\ is equivalent to exactly q + 1 distinct partially reduced 
quadratic form with the same discriminant. 



Proof. Consider the equation — {A/h)fP = 1 in part (3) of Theorem 3.3 This 
equation always has a solution, namely (a, /3) = (1, 0). Since 4//i is a non-square 
in F*, this equation is a non-degenerate conic in a and j3. Since this conic has 
at least one solution, it has q + \ distinct solutions (see Casse 9J, page 140 or 
Hirschfeld (TH], page 141). Each of these q + 1 solutions yields a binary quadratic 



form /' equivalent to /, as given in Theorem 3.3 □ 



Definition 3.5. Let f — (A, _B, C) be an imaginary or unusual partially reduced 
binary quadratic form. 

1. If\A\ < \C\, then f is called reduced. 

2. If \A\ = \C\, then f is called reduced if it is lexicographically smallest 
amongst all the partially reduced forms in its equivalence class. 



Theorem 3.6. 1. (Existence and Uniqueness) Every imaginary or unusual 
binary quadratic form over ¥q [t] is equivalent to a unique reduced binary 
quadratic form with the same discriminant. 

2. (Finiteness) There are only finitely many reduced imaginary or unusual 
binary quadratic forms with fixed discriminant D . 

Proof. Part (1) is obvious. To obtain part (2), note that if / = {A,B,C) is 
reduced, then \B\ < \D\ = \AC\. It is clear that only a finite number of triples 
{A,B,C) in Fq[t]'^ can satisfy these conditions. □ 

Following Buell fH^ , any matrix M such that f — f o M is called an automor- 
phism of /. Automorphisms of binary quadratic forms need to be considered 
in the development of the reduction theory of binary cubic forms, and can only 
occur in very specific cases if a binary quadratic form is reduced: 

Theorem 3.7. Suppose f = (A, B, C) is a partially reduced imaginary or un- 
usual binary quadratic form such that f = f o M for some M G GL2 (Fg [t] ) . If 
\A\ < \C\, thenM = ±1, or M = ±J if B = 0, where J = If\A\^ \C\, 

then M = ( _ "^.^ ^f. ) where a, f3 E F,, A = -A/h and u := dct(A/) = ±1. // 

/3 7^ 0, then B ^ 0, and a = uP{CX — A)/B G F^. Hence, such non-trivial 
transformations exist if and only if ((CA — A)/BY + \ is a square in ¥*, in 
which case there are two such non-trivial automorphisms. 
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Proof. Let A/ = ^ " ^ ^ Since / = / o Af , we obtain from equation ( 7 1 , with /' 
replaced with /, that 

iuA^ = {2Aa + -iBf - D-/^. 

If \A\ < \C\, then |Ap < Z?, so 7 = and so we obtain fro m equation ([ij that 
a GVq with = 1. Then, as in the proof of Theorem 3.3 we obtain M — ±1 
or ± J, where the latter two cases force B = 0. 

Now suppose that \A\ = \C\. Then 7 = {Auj3)/h and 5 = uahy the proof of 
Theorem |3.3[ — {^/h)l3'^ — 1 and u = ±1 is uniquely determined. If /3 ^ 0, 
then B 0. To see this, suppose /3 7^ and B = 0. Then it follows from 
equation ([T]) that 

= A{a^ - 1) + C(A2/32) A{^X/3^) + CX^P^ = A/32(CA - A), 

therefore A = AC. Since D ^¥q and / is assumed to be primitive by Proposition 
|2.1[ 4), we have the required contradiction. 

Hence B ^ 0. Again, by equation ([I]), we obtain 

= A{a^-l)-Bu\aP + C\^/3^, 
= -A\13^ - BuXaP + CX^/3^, 
= X(5{I3{CX- A) -uaB). 

Therefore, uaB — /3{CX — A) and solving for a yields a = ul3{CX — A)/B, 
as desired. To see that ((CA — A) /BY + A is a square in F*, we substitute 
a = u/3(CA - A)/B into + A/?^ = 1 to obtain 

[((CA-A)/S)2 + A]^2 ^ ^ 

So such a non-trivial automorphism exists if and only if ((CA — A)/B)'^ + X 
is a square in F*. In that case, there are two solutions for /3, and for each of 
these solutions, there is exactly one solution for a, since u is determined by 
sgn(i3) € S. This completes the proof. □ 

We note that if \A\ — |C|, the case P — is still possible and yields the same 
four matrices as for \A\ < |C|. 



4 Reduction of Binary Cubic Forms over ¥q[t] 

In this section, we describe the reduction theory for binary cubic forms. Once 
again, this theory allows us to single out a unique representative for each equiv- 
alence class of binary cubic forms, in a similar way to the reduction theory for 
binary quadratic forms. Throughout this section, —3-D is an imaginary or un- 
usual discriminant such that sgn(— 3Z?) G {1, h}. Furthermore, h and 5* are as 
in Section |3] Recall that our irreducibility assumption forces ad ^ for any 
binary cubic form / ~ (a, 6, c, d). 
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Definition 4.1. Let f = (a,b,c,d) be a binary cubic form with imaginary or 
unusual Hessian Hf = (P, Q, R). 

1. If Hf is imaginary, then f is reduced if Hf is reduced, sgn(a) £ S, and 
if Q = 0, then sgn(c?) £ S. 

2. If Hf is unusual, then f is reduced if Hf is reduced, sgn(a) G S , if 
(5 = then sgn(d) G S , and f is the lexicographically smallest among 
all equivalent binary cubic forms with the same reduced Hessian Hf that 
satisfy the previous conditions. 



The normalization on sgn(a) is needed because / and — / have the same 
Hessian. The normalization of sgn(d) in the case when Q = is required because 
in this case, Hf has automorphisms ± J with J — (^q so f — (a, b,c,d) and 
/ = (a, —b, c, —d) have the same reduced Hessian. 

In order to obtain a unique reduced binary cubic form in the case where 
the unusual Hessian is non-trivially equivalent to itself, we employ the same 
technique that was used for unusual binary quadratic forms. That is, we adopt 
the convention of choosing the binary cubic form that is the smallest in terms 
of lexicographical order as specified in Section [3] Note that it is straightforward 
to detect whether or not a Hessian has non-trivial automorphisms: by Theorem 
3.7 it simply requires checking whether or not {4R/h + P) /Q)^ — 4:/h is a square 



in F* . The lexicographical minimization condition eliminates ambiguity in case 
of non-trivial automorphisms of the Hessian. Collectively, the above conditions 
ensure uniqueness of reduced representatives: 

Theorem 4.2. Every binary cubic form f with imaginary or unusual Hessian 
is equivalent to a unique reduced binary cubic form with imaginary or unusual 
Hessian of the same discriminant. 



Proof. By Proposition 2.3 and Theorem 3.6 / is equivalent to a binary cubic 
form of the same discriminant with reduced Hessian, so we can assume without 
loss of generality that Hf is reduced. 

If sgn(a) ^ S, then replace / by -/ = / o (-/). If sgn(a) G S*, Q = 
and sgn(d) ^ S then replace / by / o J, where J = ( o -i ) • The resulting form, 
again denoted by /, has reduced Hessian and satisfies the required normalization 
conditions on a and d. 

Finally, suppose that there exists a matrix M G GL2{¥q[t]) with Hf o M = 
Hf as described in Theorem 3.7, such that f oM also satisfies the required nor- 



malization conditions. Then replace / by the lexicographically smallest among 
/ and f o M for all permissible choices of M . The resulting binary cubic form 
is reduced. Moreover, none of the above transformations change D{f) or Hf. 

To obtain uniqueness, suppose / = (a,b,c,d) and /' = [a' ,b' ,c' ,d') are 
reduced and equivalent. Then there exists a mat rix M G GL2{¥q[t]) with /' = 



foM. Set w = det(M). Then by Proposition|2.3[ Hf, ^ u^HfoM = Hfo{uM). 
Hence Hfi and Hf are equivalent, via the transformation matrix uM. By 
assumption, Hf, and Hf are reduced and hence equal by Theorem 3.6 By 



Proposition 2.3 Hf has discriminant —3D{f). 



12 



Write Hf = iP,Q,R). If |P| < \R\, then by Theorem |3.7[ uM = / or 
— / or J or — J, where the latter two cases force Q = 0. Since det{uM) = 
det(M) = u^, it follows that -u^ = 1 or —1. Thus a' = au^^ if uM = / or J, 
and a' = -a-u"^ if ^ -/ or - J. If uM = J, then -u^ = det(uM) = det(J) = 
— 1, and hence a' = au^^ = —a, which contradicts sgn(a), sgn(a') € 5*. If 
uM = —I, then = 1, and hence a' = — aw^'^ — —a, again a contradiction. If 
uM = —J, then u'^ = det(— J) = —1. In this case, Hfi = Hf implies Q = 0, 
so sgn((i), sgn((i') G S*. But then /' — u~^f o {—J) implies d' = du^^ — — d, 
a contradiction. So we must have uM — /, and hence — 1. It follows that 
/' = u^^ f o {uM) = f o I = f- This completes the proof for the case \P\ < \R\. 

If |P| = \R\, then the lexicographical minimality of / and /' forces / = 
/'. □ 

5 Coefficient Bounds for Reduced Forms 

We now present bounds on the coefhcients of a reduced binary cubic form with 
imaginary or unusual Hessian. The motivation for these bounds is to establish 
that the set of reduced binary cubic forms up to any fixed discriminant degree 



is in fact finite, which yields a result analogous to part (2) of Theorem 3.6 for 
cubic forms. This also ensures that the search procedure in Section [7] terminates 
because the search space is finite. Moreover, we seek optimal bounds on the 
coefficients of a reduced binary cubic form so that the search procedure in 
Section [T] is as efficient as possible, as smaller upper bounds give rise to shorter 
loops. 

The following equality appears in Cremona [12,, and is easily verified by 
straightforward computation. 

Lemma 5.1. Let f = (a, b, c, d) be a binary cubic form with coefficients in ¥q[t] 
with imaginary or unusual Hessian Hf — (P, Q, R). Let U ~ 2b'^ + 27a^d~9abc. 
Then 

AP^ = [/2 + 27a^D. 



Corollary 5.2. With the notation of Lemma 5.1 we have < \P\^ and 

\a^D\ < 

Proof. Since —SD is imaginary or unusual, there can be no cancellation between 
the two summands on the right-hand side of the identity of Lemma |5.1[ Hence, 
neither term can exceed \P\^ in absolute value. □ 

Theorem 5.3. Let f = {a,b,c,d) be a reduced binary cubic form over ¥q[t] 
of discriminant D with imaginary or unusual Hessian {P,Q,R). Then \a\ < 
|D|i/4^ |6| < \D\^/\ \bc\ < |L>|i/2 and \ad\ < \D\^/^. 



Proof We have \Q\ < \P\ < \D\^/'^. By Lemma 

|I?|l/2. 



5.1 



\D\ < \P\ < 
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Now an easy computation reveals U — 2bP — SaQ. This, together with 
\U\ < |P|3/2 and \aQ\ < \aP\ < |P|3/2, implies \bP\ < max{|t/|, \aQ\} < |P|3/2; 
whence 

|6| < |P|l/2 < |£,|l/4^ (14) 

For the remainder of the claim, it suffices to prove only one of \ad\ < \D\^/^ 
and \bc\ < |-D|^/^, as any one of these inequalities, together with Q — he — 9ad, 
implies the other. 



Now by ([W]), \ac\ = \b^ - P\ < \D\^^'\ If |6| < |a|, then \bc\ < \ac\ < \D\^/^. 
If \d\ < |c|, then \ad\ < \ac\ < \D\^/^. 

Finally, suppose that \a\ < \b\ and |c| < \d\. It is easy to verify that bR + 



MP - cQ = 0. Thus, \cQ\ < \dP\ implies \dP\ = \bR\. It follows from ^ that 
\bdP\ = \b^R\ < \PRl so\bd\ < \R\. This in turn implies \c\^ = \3bd + R\ < \R[ 



and hence again from (14), j&cp < \PR\ = \D\, as claimed. □ 



These bounds are indeed sharp. An example of a reduced binary quadratic 
form over Fg^ where \a\ = \b\ = \c\ = \d\ = |D|i/4 is / = {a,b,c,d) ^ (2t + 
4, 3i + 4, 3t + 3, 3t + l), where D = + + 3. 



We use the bounds of Theorem 5.3 for our tabulation algorithm. Specifically, 
we loop over all a, 6, c, d satisfying these bounds. An upper bound on |a6c(i| 
determines how often the inner most loop is entered. By Theorem |5.3[ such a 
bound is given by \D\. Hence, \D\ is an upper bound on the number of forms 
of discriminant D that the algorithm checks for membership in the Davenport- 
Heilbronn set U. 



6 The Davenport-Heilbronn Theorem 

We now briefly discuss the Davenport-Heilbronn theorem for function fields. 
The original Davenport-Heilbronn theorem [TJ] states that there exists a discriminant- 
preserving bijection from a certain set U of equivalence classes of integral binary 
cubic forms of discriminant D to the set of Q-isomorphism classes of cubic fields 
of the same discriminant D. Therefore, if one can compute the unique reduced 
representative / of any class of forms in U of discriminant D with \D\ < X, 
then this leads to a list of minimal polynomials /(x, 1) for all cubic fields of 
discriminant D with \D\ < X. 

The situation for cubic function fields is completely analogous. We now 
state the function field version of the Davenport-Heilbronn theorem, describe 
the Davenport-Heilbronn set U for function fields, and provide a fast algorithm 
for testing membership in l/l that is in fact more efficient than its counterpart 
for integral forms. 

For brevity, we let [/] denote the equivalence class of any primitive binary 
cubic form / over Vq[t\. Fix any irreducible polynomial p € Fq[i]. Analogous 
to [3 m [TO] , we define Vp to be the set of all equivalence classes [/] of binary 
cubic forms such that \ D{f). In other words, if D{f) — i^A where A is 
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square-free, then / G Vp if and only ii p \ i. Hence, / G Clp^p if and only if 
D{f) is square-free. 

Now let Up be the set of equivalence classes [/] of binary cubic forms over 
Vq[t] such that 

• either [/] e Vp, or 

• f{x,y) = X{Sx - -fyf (mod p) for some A e F,[t]/(p)*, 7, (5 e Fq[i]/(p), 
and in addition, /(7,^) ^ (mod p^). 

For brevity, we summarize the condition f{x, y) = \{5x — jy)^ (mod p{t)) 
for some 7,(5 G F,[t]/(p) and A € Fg[t]/(p)* with the notation {f,p) = (1^) as 
was done in [TU [15] . 

Finally, we set 14 = Hp^p- The set W is the set under consideration in the 
Davenport-Heilbronn theorem. The version given below appears in |25j . A more 
general version of this theorem for Dedekind domains appears in Taniguchi [31] . 

Theorem 6.1. Let q be a prime power with gcd((7, 6) = 1. Then there exists 
a discriminant-preserving bijection between ¥ q{t) -isomorphism classes of cubic 
function fields and classes of binary cubic forms over¥q[t] belonging toU. This 
bijection maps a class [f] € U to the triple of Fq(t) -isomorphic cubic fields that 
have minimal polynomial f{x, 1). 



In order to convert Theorem 6.1 into an algorithm, we require a fast method 
for testing membership in the set U. This is aided by the following efficiently 
testable conditions: 

Proposition 6.2. Let f — {a,b^c,d) be a binary cubic form over ¥q[t] with 
Hessian Hf = {P,Q,R). Let p e ¥q[t] be irreducible. Then the following hold: 

1. {f,p) = (13) if and only if p \ gcd{P,Q,R). 

2. Lf if,p) = (l^) then f e Up if and only if p^ \ D{f ). 

In addition, classes in U contain only irreducible forms; this result can be 
found for integral cubic forms in Chapter 8 of [TU], and is completely analogous 
for forms over ¥q[t\. 

Theorem 6.3. Any binary cubic form whose equivalence class belongs to U is 
irreducible. 



By Theorem 6.1 if [f] ^ U, then fix,l) is the minimal polynomial of a 
cubic function field over ¥q{t). This useful fact eliminates the necessity for a 
potentially costly irreducibility test when testing membership in U. 

Using Proposition |6.2[ we can now formulate an algorithm for testing mem- 
bership in U. This algorithm will be used in our tabulation routines for cubic 
fimction fields. The algorithm is slightly different from the one in [26]; Hes- 
sian and discriminant values are passed into the routine, rather than computed 
inside Algorithm [T] 
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Algorithm 1 Test for membership in U 



Require: A binary cubic form / — {a,b,c,d) witli coefficients in ¥q[t], its 

Hessian (P, Q, R) and its discriminant D. 
Ensure: true or false according to whether or not [/] lies in the set U. 

1: Compute in ■= gcd(P, Q, R). 

2: li £h is not square-free then 

3: return false; 

4: Compute s -~3D/{£h)'^; 

5: if gcd(s,£^f) 7^ 1 then 

6: return false; 

7: if s is square-free then 

8: return true; 

9: else 

10: return false; 



We note here that step 2 of Algorithm [T] uses part (2) of Proposition 6.2 In 
step 5, note that ifp \ in and p \ s, then \ D and using part (1) of Proposition 



6.2 yields [/] ^ Up. If / passes steps 1-6, then s is not square-free if and only 
if there exists an irreducible polynomial p G Vg[t\ with p^ \ s and hence p\ ijj- 



Using part (2) of Proposition 6.2 again, this rules out (/,_p) = (1'^). On the 
other hand, we also have p^ \ D{f), so / ^ Vp, and hence / ^ Up. Thus, if s is 
not square- free then [/] ^ Up for some p, or equivalently, [/] ^ U. Conversely, 
if s is square-free, then the primes p dividing s, and hence dividing D, occur in 
D to the first power. Thus / € Vp for all such p, proving the validity of step 7. 

Note that steps 2 and 7 of Algorithmjljrequire tests for whether a polynomial 
F G Vq[t\ is square-free. This can be accomplished very efficiently with a simple 
gcd computation, namely by checking whether gcd(F, F') — 1, where F' denotes 
the formal derivative of F with respect to t. This is in contrast to the integral 
case, where square-free testing of integers is generally difhcult; in fact, square- 
free factorization of integers is just as difficult as complete factorization. Hence, 
the membership test for U is more efficient than its counterpart for integral 
forms. One may be tempted to try a more naive approach, namely factoring 
the square part i out of the discriminant D and then testing only the resulting 
p. Even though factorization of polynomials over finite fields is much easier than 
factoring integers, this would add an unnecessary log factor to our theoretical 
run times, so we did not attempt this approach. 



7 Algorithms and Numerical Results 

We now describe the tabulation algorithm for cubic function fields corresponding 
to reduced binary cubic forms over [t] with imaginary or unusual Hessian; that 
is, cubic extensions of ¥q{t) of discriminant D where deg{D) is odd, or deg(D) 
is even and sgn(— 3D) is a non-square in F*. 

The idea of the algorithm is as follows. Input a prime power q coprime to 6, 
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a degree bound B S N, a primitive root h of F* and the set S. The algorithm 
outputs minimal polynomials for all Fg{i)-isomorphism classes of cubic exten- 
sions of Fg(t) of discriminant D such that —3D is imaginary or unusual and 
\D\ < . The algorithm searches through all coefficient 4-tuples (a, b, c, d) that 



satisfy the degree bounds of Theorem 5.3 with \D\ replaced by such that the 



form / — (a, b, c, d) satisfies the following conditions: 

1. / is reduced; 

2. / has imaginary (resp. unusual) Hessian; 

3. / belongs to an equivalence class in U; 

4. / has a discriminant D, where deg(£') < B. 



If / passes all these tests, the algorithms outputs /(x, 1) which by Theorem 6.1 
is the minimal polynomial of a triple of Fq(t)-isomorphic cubic function fields 
of discriminant D. 

The test of whether or not / is reduced in the case where Hf is unusual 
is more involved than in the imaginary case. Recall from Theorem |3.3| and 



Corollary pU] that if Hf = (F, Q,R) is the Hessian of / and both \P\ = \R\ = 
\^\D\ and sgn(P) = 1 are satisfied, then this test requires the computation 
of g + 1 partially reduced binary quadratic forms equivalent to Hf, along with 
the determination of whether Hf is the smallest lexicographically amongst the 
q + I quadratic forms. Furthermore, by Theorem 3.7 Hf may have non-trivial 



automorphisms, potentially resulting in more than one equivalent binary cubic 
form with the same Hessian. Fortunately, the proportion of forms for which 
these extra checks need to be performed is small by Lemma [7.2| below, so this 
does not significantly affect the run time of the tabulation algorithm. However, 



the fact that the bounds in Theorem 5.3 can be attained when \P\ = \R 



whereas they can not be exactly met when \P\ < \R\, make the algorithm for 



unusual Hessians slower by a factor of g, as seen in Corollary 7.3 

Checking whether a binary cubic form has a partially reduced Hessian amounts 
to checking the conditions specified in Definition |4T] Checking whether a given 
cubic form lies in U involves running Algorithm [T] A basic version of Algorithm 
2^ loops over each of the coefficients a, 5, c, d up to the bounds given in Theorem 
5.3| We omit the description here, instead giving an improved version in Algo- 
rithm [2j First, we modify the for loop on a so that the condition sgn(a) G S 
is checked first. The other improvement involves the for loop on d. Instead of 
using the for loop on d as one might do in a basic implementation, we determine 
which degree values of d lead to an odd (resp. even) degree discriminant. This 
entails computing the quantities mi through ms and the maximum of these 
values m. The values of mi through are simply the degree values of each 
term in the formula for the discriminant of a binary cubic form. If the maximum 
value m of these terms is taken on by a unique term amongst the m^ and m is 
not of the appropriate parity, the next degree value for d is considered instead of 
proceeding further with computing the Hessian and discriminant. This allows 
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us to avoid extra computations for discriminants that do not have odd (resp. 
even) degree altogether. 

For admissible values, we compute the Hessians and the discriminant as 
before, but we can compute the quantities P and ti := be, t2 '■= c? before any 
information about d is known. We also note that the seemingly redundant check 
of the conditions deg(£') < B and deg(i?) odd (resp. even) near the end of the 
algorithm may be necessary in the event of cancellation of terms when m is 
not taken on by a unique term among the rrij (since it may be possible that 
the maximum m of the degrees on the terms of the cubic discriminant satisfies 
m > B but deg(£)) < B). 



Algorithm 2 Tabulation of Imaginary (resp. Unusual) Cubic Function Fields 
(modified version) 

Require: A prime power q not divisible by 2 or 3, a primitive root /i of F,, the 

set S = {1, /i, /i^, . . . /i^'"^^/^}, and a positive integer B. 
Ensure: Minimal polynomials for all Fg(<)-isomorphism classes of cubic func- 
tion fields of discriminant D with deg(r') odd, sgn(— 3£>) € (resp. 
deg(D) even and sgn(— 3Z?) = /i), and deg(D) < B. 

1: for deg(a) < B/4 AND sgn(a) e 5 do 

2: for deg(6) < i3/4 do 

3: for deg(c) < B/(2deg(6)) do 

4: mi := 2(deg(6) + deg(c)); 

5: mi := deg(a) + 3 deg(&) 

6: for z = to B/2 - deg(a) do 

7: ma := deg(a) + deg(6) + deg(c) + i; 

8: m4 := 3deg(6) + i 

9: ms := 2(deg(a) + i) 

10: m := max{m,i , m,2, m,3, m4, mr^} 

11: if (m is not taken on by a unique term among the rrii) OR (m 

is taken on by a unique term AND m is odd (resp. even) AND 
g"' < q^) then 

12: Compute P ■.= b'^ - Sac; 

13: Compute ti := be; 

14: Compute t2 := (? 

15: for deg((i) = z do 

16: Set / := (a,6, c, d); 

17: Compute Q := ti — 9ad; 

18: Compute R:=t2 — 3bd; 

19: Compute -3£> = -3£>(/) = - iPR; 

20: if deg(D) is odd and sgn(— 3D) G {1,^} (resp. deg(D) even 

and sgn(-3£)) = h) AND deg(£>) < B AND / is reduced 
AND [/] e U then 

21: Output f{x, 1); 



Some extra routines are needed for the unusual Hessian case and are de- 
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scribed here. First, a routine called ConicSolver is used to solve the equation 
— (4//i)/3^ = 1 for a, /3 S ¥q via brute force. Each solution pair is stored in an 
array. If the Hessian Hf = (P, Q, R) of a binary cubic form / satisfies \P\ = \R 



and sgn(P) = 1, then the solutions {a,/3) are used as in Corollary 3.4 Each 



pair (a, /3) determines a matrix Ma^fi ~ 4u'p/h ua ) ' where u = ±1 is chosen so 
that either sgn[Qi) g S when Qj ^ for any of the g + 1 Hessians (P^, Qi, Ri) 
of the corresponding partially reduced q + 1 cubic forms, or sgn{di) € S when 
Qi = 0; here, = {ai,bi,Ci,di) = f o M where {P^,Qi,Ri) = Hf. ^ Hf o M, 
and M is the appropriate matrix M^^p- Each of these matrices is applied to 
each Hessian Hf under consideration, with the value u determined as described 
above. 

If the Hessian Hf is the smallest in terms of lexicographical order amongst 
the q + 1 Hessians computed, then the corresponding binary cubic form / is 
output if it lies in U. This task is accomplished via a routine called IsSmall- 
estQuad, which returns true if Hf is the smallest lexicographically amongst the 
q+1 forms equivalent to Hf, and returns false otherwise. In the event that Hf 



is the smallest and Hf has non-trivial automorphisms as given in Theorem 3.7 



then / is tested to see if it is the smallest in terms of lexicographical ordering 
amongst itself and the two binary cubic forms equivalent to / via a non-trivial 
automorphism. This task is accomplished via a routine called IsDistCuh, which 
returns true if / is the smallest lexicographically amongst these binary cubic 
forms, and returns false otherwise. If the routine IsDistCuh returns true and 
the binary cubic form / lies in lA, the minimal polynomial /(x, 1) is output. 

If the (unusual) Hessian Hf of a binary cubic form / satisfies |P| < |P|, then 
the algorithm is much simpler, since the partially reduced binary cubic form is 
in fact reduced. That is, we simply test the binary cubic form / to see if it 
lies in lA, just like in the case of imaginary Hessians. If it does, the minimal 
polynomial /(x, 1) is output. 

Tables [l] and [5] present the results of our computations using Algorithm 
[2] for cubic function fields with imaginary and unusual Hessian, respectively, 
for various values of i?, given in column 2 of each table. The third column 
gives the number of fields with discriminant D satisfying deg(Z?) < B and the 
fourth column denotes the number of reduced cubic forms whose Hessians have 
non-trivial automorphisms. The last column gives the various timings for each 
degree bound. The results in these tables extend and correct those in [26]. The 
lists of cubic function fields were computed on a multi-processor machine with 
four 2.8 GHz Pentium 4 processors running Linux with 4 GB of RAM, with a 
180 MB cache using the C-I--I- programming language coupled with the number 
theory library NTL p5] to implement our algorithm. The results for unusual 
discriminants in this table are new and did not appear in [26j . 
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Table 1: Number of cubic function fields over ¥g{t) of discriminant 

D with dcg(D) < B, dcg(D) odd 



q 


Degree bd. B 


# of fields 


# non-triv. auto. 


Total elapsed time 


5 


3 


100 


— 


0.02 seconds 




5 


2100 


— 


1.21 seconds 




7 


64580 


— 


31.66 seconds 




9 


1877260 




26 minutes, 4 sec 




11 


45627300 




9 hours, 31 min, 45 sec 


7 


3 


294 




0.17 seconds 




5 


12642 




17.48 seconds 




7 


718494 




25 minutes, 4 sec 




9 


39543210 




21 hours, 56 min, 45 sec 


11 


3 


1210 




2.61 seconds 




5 


134310 




18 minutes, 45 sec 




7 


17849810 




1 day, 2 hours, 
9 min, 28 sec 


13 


3 


2028 




7.10 sec 




5 


318396 




55 minutes, 54 sec 




7 


58239948 




6 days, 31 min, 56 sec 



Table 2: Number of cubic function fields over Vq{t) of discriminant 
D with deg(i:») < B, deg{D) even 



Q 


Degree bd. B 


# of fields 


# non-triv. auto. 


Total elapsed time 


5 


4 


280 


10 


0.88 seconds 




6 


6480 


10 


19.06 seconds 




8 


156920 


320 


12 minutes, 3 sec 




10 


4688440 


320 


7 hours, 12 min, 32 sec 




12 


117981240 


11385 


10 days, 19 hours, 
3 min, 8 sec 


7 


4 


1077 


42 


12.18 seconds 




6 


51645 


42 


17 minutes, 28 sec 




8 


2475271 


2436 


10 hours, 
38 minutes, 33 sec 




10 


138360895 


2436 


23 days, 4 hours, 
16 min, 51 sec 


11 


4 


5722 


54 


13 minutes, 5 sec 




6 


810372 


54 


21 hours, 
54 min, 59 sec 


13 


4 


11334 


304 


40 minutes, 42 sec 




6 


2240106 


304 


3 days, 14 hours 
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Timings for Algorithm [2] are compared to those of a basic version of Al- 
gorithm [2] in Table |3] for q — 5,7. Recall that this basic version simply loops 
over all all a, 6, c, d satisfying the bounds of Theorem 5.3 without a priori elim- 
inating potential unsuitable values of d as was done in Algorithm [2] The cases 
q = 11, 13 are omitted for brevity. As seen in the third column of these tables, 
the modified algorithm is a significant improvement over the basic algorithm. 
These improvements appear to get better as the degree increases, likely because 
higher degrees give fewer chances of "bad" leading term cancellations in the 
discriminant, thereby a priori eliminating more unsuitable values of d. Also, in- 
terestingly, the improvement seems more pronounced for even degrees. Again, 
this is likely due to fewer "bad" cancellations of leading terms of the discrim- 
inant in this case. For example, the terms (6c)^ or 27{ad)'^ in D can never 
individually dominate if deg(_D) is odd. 



Table 3: Basic vs. modified algorithm timings 



q 


Degree bd. 


Basic times 


Modified times 


Basic/Mod 


5 


3 


0.09 seconds 


0.02 seconds 


4.5 




4 


7.67 seconds 


0.88 seconds 


8.72 




5 


7.19 seconds 


1.21 seconds 


5.94 




6 


3 minutes, 55 sec 


19.06 seconds 


12.32 




7 


3 minutes, 20 sec 


31.66 seconds 


6.33 




8 


5 hours, 52 min, 38 sec 


12 minutes, 3 sec 


29.28 




9 


4 hours, 
7 min, 57 sec 


26 minutes, 4 sec 


9.51 




10 


5 days, 1 hour, 
38 min, 6 sec 


7 hours, 
12 min, 32 sec 


16.87 




11 


6 days, 8 hours, 
22 min, 7 sec 


9 hours, 
31 min, 45 sec 


15.99 


7 


3 


0.73 seconds 


0.17 seconds 


4.29 




4 


2 minutes, 1 sec 


12.18 seconds 


9.92 




5 


3 minutes, 56 sec 


17.48 seconds 


13.49 




6 


3 hours, 
16 min, 40 sec 


17 minutes, 28 sec 


11.26 




7 


3 hours, 12 min, 34 sec 


25 minutes, 4 sec 


7.68 




8 


10 days. 


10 hours, 


22.61 






38 min, 17 sec 


38 minutes, 33 sec 






9 


12 days, 18 hours, 
35 min, 3 sec 


21 hours, 
56 min, 45 sec 


13.97 



The value h — 2 was chosen as a primitive root for F5, Fn and F13. For F7, 
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/i = 3 was chosen. This completely determines the set S specified in Section [Sj 
These sets were {1,2}, {1,2,3}, {1,2,4,5,8} and {1,2,3,4,6,8} for Fg, F7, Fn 
and F13, respectively. 

The worst-case complexity of the algorithm, expressed in terms of the num- 
ber of field operations required as a function of the degree bound B, largely 
depends on the size of each of the coefficients a, b, c, d of a binary cubic form 
that the algorithm loops over. It follows from Theorem 5.3 that \abcd\ < \D\, 
so the number of forms that need to be checked up to an upper bound B on 
deg{D) is of order B. This idea is fully explained in the following lemmas. 

Lemma 7.1. For s G N, denote by Ts the set of binary cubic forms f = 
{a,b,c,d) over ¥q[t] such that deg{D{f)) = s, deg(a) < s/4, deg(fo) < s/4, 
deg{ad) < s/2, deg(6c) < s/2, and sgn(a) £ S. Then 



— s^q" + 0{sq'') if s is odd , 

— s^g* -|- 0{sq'') if s is even , 



as s —?> 00. 



Proof. The number of monic polynomials in Vqlt] of degree m is q"^. Hence, 
the number of pairs of monic polynomials {G,H) with deg(G') < s/4 and 
deg(G'iJ) < s/2 is 

s/4 [s/2j-m 

= E E r'^-T.r E 9" 

■m<s/i m+n<s/2 7n=0 ri=0 

- E(<zL^/^J+^-g™) = ^^.,L^/^J+0(,^/^) 

as s — > 00. 

Let / — (a, b, c, d) E J-g- Note that ad ^ 0. Thus there are (q — 1) /2 choices 
for sgn(a) and q~ \ choices for sgn(d), for a total of Ns{iq — l)'^/2 pairs (a, d). 
Similarly, the possible number of (fo, c) pairs is N^q^. Hence the total number 
of forms in J^s is 

2 32 ' ^ 

The result now simply follows from the fact that 2[s/2j = s if s is even and 
s — 1 if s is odd. □ 



Lemma 7.2. Let be as defined in Lemma \7.1\ and Eg denote the number of 
forms in such that Hf = [P, Q, R) is partially reduced and deg(P) ~ deg(i?). 
Then 

if s = (mod 4), and Eg ^ if s ^ (mod 4). 
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Proof. Let / — {a,b,c,d) G so that Hf is partially reduced with deg(P) = 
deg(i?). Then deg(Q) < deg(F) ~ deg(i?) — s/2, so s must be even. A 
straightforward calculation yields 

bQ = cP + 3aR , ScQ = 3dP + bR , (15) 

whence deg(c) — deg(6Q — 3aR) — deg(P) < s/4 and deg(d) = deg{3cQ — 
bR) — deg(P) < s/4. Hence all of a,b,c,d have degree no more than s/4. In 
fact, the formulas for Q,P,R can be seen to force deg(a) = deg(c) = s/4 or 
deg(6) — deg{d) = s/4. It follows that Eg — unless s is a multiple of 4, which 
we assume for the remainder of the proof. 

If deg(&), deg(d) < s/4 — 1, then sgn(c)^ — sgn(i?) = —h/4 forces —1 
to be a non-square in F^, or equivalently, q = —1 (mod 4). Together with 
—3 sgn(ac) = sgn(P) = 1 and sgn(a) g S, this determines sgn(a) and sgn(c) 
uniquely. Thus, this accounts for at most (7" forms if g = —1 (mod 4) and no 
forms a q = 1 (mod 4). 

If deg(a), deg(c) < s/4 — 1, then deg(&) = deg(d) — s/4. In this case, 
P ^ b^ — 3ac forces sgn(6) — ±1, so there are 2q''/^ possibilities for b. Then 
sgn(i?) = —3 sgn{bd) determines sgn(d), so there are q^^^ choices for c and d 
each. The number of permissable a is at most 



m=Q 



So this case produces no more that forms. 

Finally, suppose that a,b,c,d all have degree s/4. By (15), the leading co- 
efficients of a and b determine those of c and d uniquely. Specifically, sgn(c) = 
3/1 sgn(a)/4, and substituting this into P— 6^— 3ac yields sgn{b)^—h{3 sgn(a)/2)^ 
1. The equation — hu^ = 1 represents a non-degenerate conic which has q + 1 
solutions {u, v) over Fg. However, since a 7^ 0, the two solutions (u, v) — (0, ±1) 
are invalid. In addition, if q = —1 (mod 4), then v — produces two invalid 
solutions. In fact, sgn(a) G S eliminates half of the remaining solutions, allow- 
ing {q — 3)/2 choices for (sgn(a), sgn(&)) when q = —1 (mod 4), and {q — l)/2 
choices when q = 1 (mod 4). So the number of possibilities for / is (q — 3)q'^ /2 
if g = -1 (mod 4) and {q - l)gV2 if g = 1 (mod 4). 

In all scenarios, the number of possibilities for / adds up to {q + l)q^ /2 as 
claimed. □ 

Corollary 7.3. Assuming standard polynomial arithmetic in Vq[t], Algorithm 
2 requires 0{B'^q^) — 0(g^+'^) operations in ¥q as B ^ 00. The O-constant is 
cubic in q when B is odd and quartic in q when B is even. 

Proof. Algorithm [2] loops exactly over the forms in J^g for s < B, with Tg 
as in Lemma |7.1[ For each such form /, the entire collection of polynomial 
computations in Algorithm [2] including those of Algorithm [l] requires at most 
Ks^ field operations for some constant K that is independent of B and q. This 
holds because all polynomials under consideration have degree bounded by s. 
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The overall complexity of Algorithm [2] is dominated by step 20. The sign and 
degree checks in that step, as well as the test for partial reducedness of Hf = 
{P,Q,R), require negligible time. When \P\ — \R\, the q+l partially reduced 
forms equivalent to Hj must be computed and compared to Hf. Once the 
reduced Hessian is found, one needs to check if it has non-trivial automorphisms 
according to Theorem |3.7| If yes, compute the two corresponding equivalent 
binary cubic forms with the same Hessian and identity the lexicographically 
smallest among these three cubic forms. It follows that the overall complexity 
of Algorithm 2 is certainly bounded above by 

B 

s=l 

By Lemma pT^ 2{q+l)Es — 0{q^), so this contribution is neghgible compared to 
the error term in the bound on #J^s given in Lemma [7J] Hence, the asymptotic 
complexity of Algorithm 2 is bounded above by 



B 

E 

s— 3 s— 3 



{{Css'q' + 0{sq')) ■ Ks') < J2 Csq' ■ KB* + 0{B^q^) 



where Cg — 'J''^/32 if s is odd and Cg = 'i''*/32 is s is even. 
Suppose first that B is odd. Then 

B (B-l)/2 3 (S-l)/2 4 „5 

s=3 i=l i=2 ' 

Similarly, if B is even, then 

B B/2-1 3 B/2 4 

i5=3 i=l i=2 

< q- + qB = + 1) . 

32(^2-1)'^ 32((72 -i)^ 32(g2-i)'' 
Hence, the overall run time is 0{B'^q^) where the O-constant is as claimed. □ 



The bounds in Corollary 1 7 . 3 1 appear to be reasonably sharp. For any fixed q, 
the time required to run Algorithm [2] on the discriminant degree bound B + 2 
(of the same parity) should be larger by a factor of q^, compared to the time 
required when using the bound B. If _B is odd, then the computation time 
of Algorithm [2] using the (even) discriminant degree bound B + 1 should be 
larger by a factor of (1 + B~^)'^{q^ + l)/2, compared to the time required when 
running the algorithm with the bound B. Similarly, going from an even bound 
B to -B + 1 should increase the run time by a factor of (1 + B"!)" • 2q^/{q^ + 1). 
Our computations times in Tables 1 and 2 largely bear this out. 

If we write q^ = X, i.e. \D\ < X, then the complexity of Algorithm 2 is 
0(X^+'^) as X oo, which is completely analogous to the run time of Belabas' 



24 



Log(Time) vs. Degree Values in F5 
Line of best fit (imaginary) f(x)= 1.1 1451*x - 5.62294 ■ 
Line of best fit (unusual) g(x) = 1.08620*x - 4.55001 



log_5(X) 



Figure 1: F5 Tabulation Timings for Cubic Function Fields: logs (sec) versus 
degree 



algorithm [3] for tabulating cubic number fields of absolute discriminant up to 
X. 

Corollary |7.3| states that our algorithm in the imaginary and unusual cases 
should be roughly linear in X if 5 is small. To see that this is the case in 
practice, we plotted the various values of log^(sec) versus \og^X (i.e. degree) 
for q = 5 and q — 7, where sec denotes the time (in seconds) taken to tabulate 
all cubic function fields whose discriminant has absolute value at most X for 
various values of log^ X. The line of best fit for the data in the imaginary and 
unusual Hessian cases is also given in each figure. As seen in Figures [l] and [2] 
the running times (in seconds) of Algorithm [2] are approximately linear in X for 
both imaginary and unusual Hessian, as expected. 

The algorithm presented in this section has some of the same advantages 
as Belabas' algorithm [3] over earlier field tabulation algorithms (see Cohen 
[TU], Chapter 9). In particular, by Theorem 6.3 there is no need to check for 
irreducibility over ¥q(t) of binary cubic forms lying in U, no need to factor 
the discriminant, and no need to keep all fields found so far in memory. Our 
algorithm has the additional advantage that there is no overhead computation 
needed for using a sieve to compute numbers that are not square-free, by the 
remarks following Algorithm [l] 

The number of binary cubic forms with Hessian having non-trivial automor- 
phisms was rare. The percentage of such cubic fields with discriminant degree at 
most B having non-trivial automorphisms appears to be tending towards zero 
as B — >■ 00. We conjecture that this rare behavior persists for higher degrees. 
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Log(Time) vs. Degree Values in F7 
Line of best fit (imaginary) f(x) = 1.12036*x - 4.19351 ■ 
Line ofbest fit (unusual) g(x)= 1.01827*x - 2.69310 ■ 



log_7(X) 



Figure 2: F7 Tabulation Timings for Cubic Function Fields: log7(sec) versus 
degree 



but analyzing the behavior of such fields and their distribution remains an open 
problem. 

8 Conclusion and Open Problems 

The main results of this paper are the development of the reduction theory of 
binary cubic forms with coefficients in ¥q[t\, and its use in conjunction with 
the Davenport-Heilbronn theorem to obtain an algorithm for tabulating cubic 
function fields. The tabulation algorithm checks 0{q^^^) reduced forms in order 
to tabulate all cubic function fields with imaginary or unusual Hessian whose 
discriminant satisfies deg{D) < B, which is in line with Belabas' result [5] for 
number fields when q is small. 

The reduction theory developed here is applicable to cubic forms (and hence 
function fields) of discriminant D when —3D is imaginary or unusual. It is 
unclear which suitable quadratic form should be associated to a cubic form 
of discriminant D when —3D is real. Neither the Mathews [23j, Berwick and 
Mathews [5] , nor the Julia approach [T^] appear to be applicable in general here; 
even if there are certain cases where they might lead to a unique representative 
in each equivalence class of cubic forms, it is unclear how to derive upper bounds 
on the coefficients of such a form, due to the non-Archimedian nature of the 
absolute value on [t] . 

One possible way to overcome this obstacle is to change the question some- 
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what. Instead of considering cubic extensions F/¥q{t) of discriminant D up 
to some degree bound, we consider such extensions whose ramification divi- 
sor (or different) D has a norm which satisfies the degree bound. Here, D 
incorporates the information on all the ramified places, including the infinite 
ones, while only the finite places are contained in the discriminant D. We have 
deg(A/'(S))) = deg(D) + ep where ep (0 < e_F < 2) is given by the ramified 
infinite places of F and can be computed from the signature at infinity of F 
(see [20l [211 [28]). Here, one needs to understand the relationship between the 
ramification divisor of a cubic form and that of a suitable associated quadratic 
form. A more detailed exploration of this approach is the subject of future work. 

From our tabulation output, it appears that the number of cubic extensions 
over Fg(t) with odd discriminant degree is always divisible by q{q — 1). For 
imaginary Hessians, the divisibility by q is easily explained: every one of the q 
translates t ^ t + u with u S F* keeps a form reduced since it does not change 
any degrees or signs. The resulting form is different unless a, 6, c, d are all 
polynomials in f — t, which is impossible from the degree bounds in Theorem 



5.3 for reasonably sized D: unless deg(£') > Aq, a and b must be constant, 
which forces deg(c) = deg(d) > q. Then deg(£') = 3deg(c) > 3g, which is still 
very large. For unusual Hessians, the above argument does not apply is as the 
lexicographical ordering would not be preserved under these translates. It would 
be interesting to be able to prove if this type of divisibility phenomenon always 
occurs, or at the very least, prove specific formulas for fixed even discriminant 
degree and q values. We discuss this in an upcoming paper fI7\ . 

An explicit comparison to the Datskovsky- Wright asymptotics on cubic func- 
tion fields 13] was not completed here, since we did not consider the case where 
—3D is a real discriminant. Furthermore, the asymptotics are not given for 
each possible signature for —3D. Other asymptotics on function fields of ar- 
bitrary degree which take into account the Galois group include Ellenberg and 
Venkatesh |16j . Density results for number fields can be found in [G] I17j , among 
others. 

Constructing tables of number fields has been done for cubic, quartic and 
other higher degree extensions (see Cohen [TO])- To the knowledge of the au- 
thors, the problem of tabulation of function fields has not been widely explored. 
The generalization of existing algorithms used for tabulating number fields to 
the function field setting is also the subject of future work. 

Belabas modified his tabulation algorithm to compute 3-ranks of quadratic 
number fields W . This has also been generalized to quadratic function fields in 
pS] . and is the subject of a future paper |27j . 
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